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Symmanzik equation and in the limit the field goes to zero 
recover the known results. The effect of a finite temperature 
on the system is discussed at the end of the paper. We sug- 
gest to use the measurement of x(h) to determine the ex- 
change interaction J of ID Heisenberg antiferromagnets. 

We start bosonization by using the Jordan-Wigner trans- 
formation which maps the spins to fermions. 33 The Heisen- 
berg model H—J^ljSj- Sj +i + gfi B H^,jSj z becomes 



#=-^2 (c] +1 Cj+c}Cj +1 )+jJ, 



i\ ^ / i 

+ 8/^bH2j \rij- 



(1) 



The anticommuting operators Cj are for the pseudofermions 
and rij the number operators, rij=CjCj. In the absence of 
the field H, the free part of the above Hamiltonian describes 
a fermion system at half filling with k F — irl(2a). From the 
spin viewpoint, this means that the ground-state magnetiza- 
tion per site M z —Q. In the presence of the density-density 
interaction, the picture of Fermi liquid is not quite applicable 
in ID and should be replaced by that of a Tomonaga- 
Luttinger liquid. 34 But the quantity k F still has its signifi- 
cance not only as the starting point of bosonization, but also 
as a measurable quantity. It reflects the spatial periodicity of 
the staggered SDW's which enters the dynamical suscepti- 
bility as the locations of the ± 2k F peaks when the frequency 



is zero. ' It is also related to the change of particle number, 
and hence the magnetization, by 



M z (H)= 8J ^Vk F {H)-k F l 



(2) 



A finite field H shifts the fermi surface away from half filling 
and the magnetization is then finite. The vacuum of the sys- 
tem is accordingly changed such that the normal order of the 
operators rij is now defined as :nj--=n—k F {H)lTT. With 
the above changes, we can use the standard operator 
transformation to bosonize the fermionic fields: 
Cj^ip(x)^e ik F x t// + (x) + e~ ik F x t//_(x), with tp+{x) 
= (V^2Tra)e i ^' ,, + (x) and <M*) = ( H^Tra)e "'^-M 
In the expressions we have taken the continuum limit and 
represented the coordinate of the j'th site to be "x." The 
operators 4> ± are the chiral bosonic fields. For details of the 
procedure see Ref. 19. The bosonized Hamiltonian is 



»=l\ 



dx 



K o (d x 0) 2 +^(d x (f>) 2 



+ g j dxcos[qx+4l(m(j>(x)] (3) 

with g = 2Ja/(2ira) 2 as the coupling constant, v=Ja{\ 
+ l/7r)Vl-<i 2 , K =^l(l-d)/(l+d), and d=(2/l 
+ 7r)sin(&fa). The momentum q in the Umklapp term repre- 
sents the deviation of the fermi surface, q = 4k F — 4k F . The 
SDW speed v and the boson stiffness K are in agreement 
with those given in Ref. 33 when H—0. In a weak magnetic 
field their dependences on H are found to be second order. 
When the field H—0, and so q — 0, the Hamiltonian Eq. (3) 
has been studied using the RG method. 35 The initial values 



of the scaling variables, x = 4K — 2 and g , were found to 
reside on the critical line and scale to the fixed point, 
(x*,g*) — (0,0). The initial value of x is shifted by the field 
such that 



x,5(7/ ) -v,]- , 



(4) 



where h\ is a positive constant. The Hamiltonian Eq. (3) is 
defined at the microscopic length scale, the lattice spacing 
"a." The field 6 is dual to the bosonic field 4>. These fields 
obey periodic boundary conditions 6(x + L)= 6{x),4>{x 
+L) — (f)(x), where the length of the system L^c°. The 
coupling constant g is defined on the lattice with momen- 
tum cutoff A = 2-7r/a. When the cutoff A is reduced to 
Alb(b>\) and an integration in the momentum shell 
[A/Zj,A] is performed one generates new terms of the form 
\d x [<p(x) + qx]} 2n (n^2) and cos{2m(4l6Trc/>(x) + qx]}(m 
2= 1 ). As a result we find that the renormalized Hamiltonian 
at the scale b > 1 takes the form 



v(b) f 1 
H(b)= — \ dx K(b)(d x 6) 2 + -— {d^) 2 



K(b) 



+ 8o(b) j dx cos[(qx)z(b) + Vl6-7r0(x)] 
+ 2 (\ 2( „ +1 )(fc)(^) 2( " +1) 

n = 1 

+ g 2n (b)cos{2n[{qx)z{n)+ 4^yn<t>{x)']}) 

GO 

+ 2 X 2m ,2n(b)M) 2n cos{2m[(qx)z(b) 

m,n — l 

+ Vl67r^(x)]}. (5) 

The Hamiltonian in Eq. (5) is characterized by: z{b—l) 
= 1 and A 2 («+i)(b=l) = g2«(b=i) = £2m,2«(>=;L) = (). We 
introduce the new variables x(b) — 4K(b) — 2 and y(b) 
'xgo(b) and find the following RG equations for the cou- 
pling constants: 



dx 
dlnb 



= -y 



dy 

d)nb 



= -xy, 



dz 
d]nb 



= -cy 2 , (6) 



where c — ^21 it. The proportionality constant in defining 
y(b), which is nonuniversal, is chosen in such a way that the 
coefficient on the right side of the first equation is of mag- 
nitude "1." The rest of the parameters have negative dimen- 
sions and are therefore irrelevant, d]n\ 2 ( n+ i)/dlnb<0, 
d\ng 2n ld\nb<Q, and d\n\ 2m2n ld\nb<Q. In addition to 
Eqs. (6) we have two more relations: v K — v(b)K(b) and 
the scaling equation for the momentum q {^H to the linear 
order in the magnetic field), q(b) — bqz(b). We investigate 
the flow in the gapless phase x 2 (b)—y 2 (b)^0 for weak 
magnetic fields h< 1 . From the scaling Eqs. (6) we observe 
that q(b) grows with the scaling. At the scale b = b such 
that 



q(b ) = b qz(b ) = Tr/a, 



(7) 
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the Umklapp term becomes irrelevant. This corresponds 
to q(b ) = 4[k F (H)-Tr/(2a)] = TT/a, or k F (H) = ?>irl 
(4a). At this point, the Umklapp interaction is 

l) ; 'cos[-\/l67r0(x J -)] oc g(feo)/ dx sin[\Jl6ir(/>(x)]d x (t>. 
For scales b^b we obtain the following model (s^ 1): 



= W 



dx 



K(s)(d x e) 2 +^—(d x <p) 2 

K(s) 



+ G(s) \ dxsin[yjl6iT<f>(x)]d x <f> 



+ G 2 (s) dxcos[2^ 



+ \ 2;2 (s) dxcos[2^16TT(f>(x)](d x 4>y + 



(8) 



where G(s = 1) = g(b )^fl(m, K(s=l) = K(b ), G 2 (s=l) 
= g 2 (b ), and \ 2;2 (s= 1) = \ 2;2 (Z> ). From Eq. (8) we ob- 
serve that the new coupling constants G, G 2 , and X 22 are 
strongly irrelevant, which result in the absence of scaling of 
K(s), dK(s)/d]ns=0. Hence we identify a continuum set of 
fixed points K(s) s ^ x ^K(b ), controlled by the value of b 
and the initial value of the magnetic field h. As a result we 
obtain the longitudinal magnetic susceptibility x z W 
= 2( gf i B ) 2 K(b )/(j7T 2 ) where tf(6 ) = £[l + £x(6 )], or 



Jtt z 



1 

1+2^0) 



(9) 



In order to clarify the singular behavior in Eq. (9) we restrict 
ourselves to the critical line x(b=l)=x and y(b—l) 
= v , such that x 2 — y 2 = 0. This is valid by neglecting the 
second-order field dependence in Eq. (4). Integrating Eqs. 
(6), we obtain 



x(b ) = 



x 



l+x ln(fc ) - 
From the conditions q(b ) — irla and q^h we find 



Hb ) 



l- 



CXq 



l+x ln(fc ) 



(10) 



(ID 



Equation (10) can be solved in two critical regions: (i) h 
—s-0 such that x ln(fo )>l. From Eq. (11) we obtain ]nb 
=ccb+ln(l//j)=ln(/io//i). Hence x(Z> ) = l/lnZ7 =l/ln(/z /%)- 
(ii) The nonuniversal region, x ln(b )<l. From Eq. (11) we 
obtain lnb — (1— cx^) _1 ln(l//z), giving rise to x(b )— x {l 
+ [x /(l -cxo)]lnft}— xq/j'' where v =x /(l -cxq)>0 is 
a nonuniversal exponent. We introduce the crossover mag- 
netic field I=exp(-l/x ) (which depends on x ), and ob- 
tain the susceptibility 



Xz( h )= T 2 



1 

1 + - 



1 



2 \n(h /h) 
1 + ^"° 



h<h 
h>h, 



(12) 



1.20 



1.15 



1.10- 



1.05- 



1.00 -t 




FIG. 1. The plot of m(h)/h vs h. The solid line is drawn using 
Eq. (13) at h =6.8. The crosses are the numerical data taken from 
Ref. 25. 

where /z =exp(cx ). We can now compare our susceptibility 
with the numerical data given in Ref. 25. We consider the 
limit h^0 where the lowest-order logarithm contribution is 
the major correction. Integrating the above equation in the 
region h<h, we obtain the dimensionless magnetization 



m(h) = 



Jtt 1 



h 

dhx z (h) = h+ —Ei 



My) I (13) 



where £i(x) is the exponential-integral function defined as 
E { (x) = f^e~ x 'dt/t for x>0. We find the nonuniversal con- 
stant h — 6.8 (see Fig. 1). From the figure we see that the 
expression gives the asymptotic behavior of mlh as the filed 
h goes to zero. To get more precise value of the constant h , 
magnetization data at fields closer to zero are needed. 

We are now in the position to calculate the spin-spin cor- 
relations in the magnetic field. We concentrate on two limits. 
The first is the infrared hmit where the length scale R>b 
~h lh. In this hmit the problem is controlled by the free 
boson spin liquid with the boson stiffness given by the fixed 
point value K(bo). As a result the transversal spin correla- 
tion is determined by the correlation of the boson field 6 with 
the exponent l/[2K(b )]. The longitudinal part is controlled 
by the bosonic field </> with the exponent 2K(b ). So we 
have 



(S + (R)S.(0)) R>b 



cos[2k F (h)aR] 

" /jl-l/[2 Hh /h)] 



const, 



([S z (R)-M z ][S z (0)-M z ]) 



R>b n 



cos[2k F (h)aR] 

' ^l + l/[21n(A /ft)] 



const, 



(14) 



where the Fermi momentum k F (h) is related to the magne- 
tization by Eq. (2). The second is the nonuniversal hmit R 
<b . This region is in principle not accessible by the RG 
since the irrelevant coupling constants, which are not in- 
volved in the RG Eqs. (6), contribute. We will treat them 
perturbatively by keeping only the most important contribu- 
tions, those from the Umklapp interaction in Eq. (3), in the 
following calculation. We also make use of the fact that un- 
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der this limit, the magnetic field can be ignored in the Um- 
klapp term in Eq. (3). The only effect of the magnetic field is 
given by the shift of x to x (h) given by Eq. (4). Hence we 
can use the regular sine-Gordon scaling equations [the first 
two in Eqs. (6)] on the invariant fine x 2 (b) — y 2 (b) 
= (h/h l ) 2 , for b<b . The RG flow is found to be 

h I h \ hlh x 

(15) 

Using the Callan-Symmanzik equation for the spin-spin cor- 
relations, we have 



{S+(R)S-(0))^cos[2k F (h)aR] 

'Mb 

~b 



Xexp 



r 

Jo 



([S z (R)-M z ][S z (0)-M z ]) 



— cos[ 2 k F ( h ) aR ] exp 



Rdb 



Mb) 



l-y(b) 



const, 



+ ^x{b) 



const. 



(16) 



For details of the derivation, see Ref. 28. Substitute the RG 
flow, Eqs. (15), into the above expressions, and we get 



{S + (R)S.(0)) R<bo 



cos[2k F (h)aR] 
R 



sinh -r- \aR 
\h x 



1/2 



const, 



{[S z (R)-M z ][S z (0)-M z ]) R<bo 



cos[2k F (h)aR] tanh[(/i/2/i 1 )ln/?]/(/i/2/i 1 ) 



R 



[sinh[(^ 1 )lni?]/(^ 1 )] 



1/2 



const. 



(17) 

When the magnetic field h vanishes, the anisotropy between 
the transversal and longitudinal correlations disappears, re- 
producing the critical sine-Gordon correlation 
(-l) R \]]nR/R. The result can also be obtained directly by 
calculating Eqs. (16) on the separatrix line x 2 (b)—y 2 (b) 
= 0. 

Finally we want to make some remarks for the cases when 
the temperature is finite. Any practical measurement of the 
magnetic responses can only be done at finite temperatures. 
So far we have discussed only the field dependences of the 
physical quantities at T=0. The results are still correct as 
long as the temperature is low enough such that h /h 
<T /T where T is a constant. On the other hand, in the 
weak-field limit T /T<h /h, the temperature plays the role 
to cutoff. We need then to replace b ~h /h everywhere in 



the scaling we discussed above by T /T. The replacement 
results, in the limit h—>0, ir 2 Jx z (T)l(giJL B ) 2 =\ 
+ 1/[2 ln(r /D] which was found in Ref. 14. The RG flow in 
the weak-field limit is not critical. Only at length scales 
shorter than the correlation length determined by the inverse 
temperature does one observe critical-like behavior, while in 
a finite magnetic field the flow actually results in a new fixed 
point with an infinite correlation length. 

One possible application of the results presented in the 
paper is that the asymptotic expression of the susceptibility, 
Eq. (12) in the weak-field region, can be used to measure the 
exchange energy J of ID Heisenberg antiferromagnets. The 
system can be kept at a low but finite temperature which is 
above the 3D ordering temperature. When an increasing 
magnetic field is applied, we will first see the region domi- 
nated by the temperature cutoff where x(h) changes slowly. 
After the field becomes strong enough to play the role of 
cutoff such that h Ih < T IT, we are in the regime where the 
asymptotic expression of xW describes the magnetic re- 
sponse as long as h is small. Since xW will increase much 
faster than that in the first region, we can easily tell the 
crossover between the two regimes. For an exchange energy 
7~2000 K, a field as strong as H— 10 T corresponds to a 
dimensionless field /z~7X10 -3 , which is deep enough in 
the asymptotic region such that our expression of ^(/z) is a 
good approximation. 

In conclusion, we analyzed the response of an AF Heisen- 
berg chain in external magnetic field. We got the perturba- 
tive RG flow equations. The RG flow was found to start near 
the critical fine of the h = case and later, at the length scale 
h /h, cross over to a region with the fixed point described 
by a free boson system. We computed the uniform suscepti- 
bility x z (h) m & compared it with the numerical result. The 
infinite slope of XzW as the field goes to zero was found to 
originate from a logarithm contribution: XzW^i^ 
+ 1/[2 ln(ft //j)]}. We calculated the spin-spin correlations. 
Our results are compared with those at finite temperatures in 
the limit of zero magnetic field. We suggested to determine 
the exchange energy J by the measurement of the suscepti- 
bility x z (h). 

Note added. After our paper was submitted, we learned 
about a paper 36 by Affleck and Oshikawa which addresses 
similar problems, in particular the susceptibility of the 
Heisenberg model in an external magnetic field. Although 
the method used by us is different from that by the authors, 
the results are similar. The main difference is in the proof of 
the existence of the set of fixed points by us. 

The authors would like to thank Professor Joseph L. Bir- 
man and Professor B. Sakita for interesting discussions and 
for reading the manuscript. P.S. wishes to thank Professor V. 
P. Nair and Professor Stu. Samuel for many useful conver- 
sations. 
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